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Abstract: Cancer remains a significant burden for HIV-infected individuals. In this paper, an ODE 
model and an DDE model about HIV-1 dynamics incorporating the AIDS-related cancer 
cells in vitro a studied. For each model, we discuss the existence, the stability properties and 
the biological meanings of these steady states. Also we find conditions for Hopf bifurcation 
of the positive steady state, leading to periodic solutions, sequences of period doubling 
bifurcations and appearance of chaos. Further, chaos and periodic behavior alternate. 
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1 Introduction 


Cancer remains a significant burden for HIV-1 infected individuals. Most of these are 
virus-associated cancers. How can the combination of immunosuppression and activation of 
inflammation promote cancer development? Our purpose in this paper is to try to give a 
glancing analysis using two simple dynamical models. 

According to the literaturel"], the cell-to-cell mechanism of HIV-1 transfer has been esti- 
mated to be much more important than infection by free virus particles in lymph tissues, where 
98% of the CD4+ T cells in vivo are found. Therefore, an ODE model concerning the cell-to-cell 
spread of HIV-1 is relevant. 

There is a lag between viral infection and viral production: for a normal activated infected 
CD4+ T cells, it is estimated that the first viral release occurs about 24 hours after the initial 
infection. While for some quiescent infected lymphocyte cells, the “incubation” phase can last 
about 60 days!?]. So the second aim of the present paper is to focus the attention on the 
“Incubation” function during the HIV-1 infection and cancer developing. 

The basic starting point of this model has three parts. First, the cancer is a single clone 
mechanism. Second, the cancer cells have some special genes and so they proliferate in a 
special way which is different from normal cells. Third, considering the key position of CD4+ 
lymphocyte, we use them to represent the immune system in our model. 
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2 The ODE model 


We use C(t), T(t) and I(t) to represent the concentrations of cancer cells, healthy cells and 
infected cells respectively. We develop an ODE model as follows 


dc = C(e) |r (1 - SOH ) -kT()], 


dt 
T STO [ra( 1- SO ) — pk, Cle) = klt) ], (1) 
dI 


We use parameter rı to represent its uncontrolled proliferation rate, kı as the immune 
system’s killing rate of cancer cells; m as the effective carrying capacity of the system; r2 as 
the intrinsic growth rate of healthy cells and p as the losing rate of the immune cells because 
of its killing the cancer cells. k2 is the infection rate coefficient that accounts for the overall 
effects of HIV-1 reproduction. jz; represents the whole immune system’s killing effect on the 
infected cells. All the parameter values above are non-negative. We assume initial conditions 
of C(0) = Co, T(0) = To, I(0) = Ío. 

The possible steady states of system (1) are as follows: the trivial steady state, Eo = (0, 0, 0); 
the cancer steady state, E, = (m,0,0); the healthy steady state, E2 = (0, m, 0); the cancer- 
healthy steady state 
MT2 mMTıp ); 


Ea Os t +p(rı + mki)’ r2 +p(rı + mkı)’ 0 


the HIV-healthy steady state 


? 


Hr T2(mk2 — pr) ) 


m=0f d= (0.2, eee a 


and the cancer-HIV-healthy steady state, E* = (C*,T*, I*), in which 


mr k3 — (riko + kir2 + mky ke) pr T* = pI 


C*= 
ry ko(ke — pki) k2 
ra ki[pmrik2 — (r2 + pri + pmkı)uz] 


rıka(pkı — k2) 


About the existence and stability of these steady states we have the following results. 
Theorem 1 Let 
k k -1 -1 
Ro= 2. (sa 2(Ro ) pes i; (2) 
HI ky mk2z+r2 r2 + mpkı 
1) Æo and £3 always exist and are unstable; Æ; always exists and is locally stable. 
2) When Ro < 1, Ez is locally stable, E4 and E* do not exist. 
3) When Ro > 1, Eo is unstable and E; exists. When R; < 1, Ex is locally stable; E* 
does not exist or if it exists is unstable. When R, > 1, E4 is unstable. When Rọ < 1, E* exists; 
otherwise, E* does not exist. 
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Theorem 2 1) When aiaz — a3 > 0 holds, here 


TL ko 


Q1 = io + age. a3 = (ke — pki) —C*T"I*, 
m m m 


a2 = ka( ow kp JTI = in (p+ at) + pky err", 
m m m 
then the infected steady state E* of system (1) is locally stable. 

2) Otherwise, choose parameter rı as the crucial parameter. Then a Hopf bifurcation will 
occur when rı passes through the critical value rf, where rf satisfies: a pair of purely imaginary 
roots exists for the characteristic equation of E* when rı = rf; and the transversality condition 
q4 Reda,3|r:=re # 0 holds. 

We omit the proof here. In biology, when the uncontrolled proliferate rate of the cancer cells 
rı is small, the cancer situation can not be built and only HIV can survive in the individual. 
Otherwise, the cancer structure will develop soon, possible denoting the AIDS phase in the 


infected individual. 


3 The DDE model 


According to the immune response caused by HIV-1 in tissue culture, we employ a fixed 
delay 7]. We propose a DDE model as follows 


£69 op [n (1- arosa) nro], 


dt 
oe = T(t) [r(1 - gormet ) -— pkıC(t) — ka1(t) |, (3) 
EN = kT (t — TM (t — rT) — pr (t). 


We assume constant initial conditions: C (0) = Co, T(@) = To, I(0) = Io, for any @ € [—7, 0]. 
The existence of steady states of system (3) is the same as that for system (3). The stabilities 
of these steady states will be shown as follows. 
Theorem 3 Let Ro, Rı and Rz as that in equations (2). Then for the delay system (3), 
we have Eo and £3 is unstable for all r > 0; E; is absolutely stable for all r > 0. 
When Ro < 1, E is absolutely stable; otherwise, E2 is unstable for all r > 0 and E; exists. 
(a) When R; <1 
(i) If 3u7 > mkz holds, then E4 of system (3) is asymptotically stable for all 7 > 0; 
(ii) If 3ur < mkz holds, then E4 is asymptotically stable for r < ĉo, and unstable 
when r > ĉo, where 
ee (Ag — A1A3)@§ — A244 
= Do arccos ( area e | aaa + AB J 
When 7 = fo, a Hopf bifurcation occurs; that is, a family of periodic solutions bifurcates from 
E, as T passes through the critical value 7. 
(b) When Rı > 1, E4 is unstable for all r > 0. When R2 < 1, E* exists; suppose 
a1a2 — a3 > 0 holds, then the infected stable state E* of system (3) is asymptotically stable 
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when 7 < 79, and unstable when 7 > To, where 


1 
To = — - arccos 


( bswg + [bi (be = baw?) = bobs|we am: bz (be = baw) ) 
wo 


(bg = baw?)? + bwg 
When T = 70, a Hopf bifurcation occurs; that is, a family of periodic solutions bifurcates 


from E* as 7 passes through the critical value 7. 


4 Numerical simulation and discussion 


We mainly discuss the mathematical behavior with changing rı and kı. We choose other 
parameter values as r2 = 0.03, ko = 5e7*, m = 1500, ur = 0.3, p = 0.1, then we can obtain 
the Hopf bifurcations regions for system (1) and (3) shown in Figure 1 (the above two figures). 
Simulations for system (1) show that, the orbitally asymptotically stable periodic solution and 
chaos appear alternately with rı increasing (the bottom two figures). 
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Figure 1: The Hopf bifurcation regions of the ODE model and the DDE model, respectively 
(the above two figures). Period 1 orbit and chaos appear for rı = rg +8 x 10-5 and rı = 
rj + 0.0235 respectively for the ODE model (the botton two figures). 
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Considering the HIV-1 infection in vivo, i.e., considering HIV virus inside the model, we 
still can find similar results as the two models before. In fact the existence of oscillations and 
chaos in a biological model is not novel. This phenomenon has been found in clinics!4!. Also, 
the research on the throb of the human heart and the galvanic activity of the human brain 
found chaos phenomena‘). 
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